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We provide a full quantum mechanical analysis of a weak energy measurement of a driven me- 
chanical resonator. We demonstrate that measurements too weak to resolve individual mechanical 
Fock states can nonetheless be used to unambiguously detect the non-classical energy fluctuations 
of the driven mechanical resonator, i.e. "phonon shot noise". We also show that the third moment 
of the oscillator's energy fluctuations provides a far more sensitive probe of quantum effects than 
the second moment, and that measuring the third moment via the phase shift of light in an optome- 
chanical setup directly yields the type of operator ordering postulated in the theory of full-counting 
statistics. 
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Introduction. — There has been considerable interest 
recently in preparing and detecting quantum mechanical 
behaviour in mechanical resonators. This general goal 
has been pursued using both optomechanical systems 
(where the mechanics is coupled to an optical cavity) 
and electromechanical systems (where the mechanics is 
coupled to an electrical circuit). Not only do such studies 
attempt to test quantum mechanics in a new regime of 
large mass, they also have the potential of furthering our 
understanding of quantum dissipative processes and the 
boundary between classical and quantum physics. Per- 
haps the simplest example of truly quantum behaviour is 
the energy quantization expected of a quantum oscillator. 
Detecting this quantization directly by, e.g., observing 
quantum jumps between different Fock states [1], is ex- 
tremely challenging [SHI]. One requires a detector which 
couples directly to the energy of the mechanical oscilla- 
tor, not to its position x as is more common; further, this 
coupling must be strong enough to resolve a mechanical 
energy quantum within the lifetime of a Fock state. Re- 
cent optomechanical experiments [4] have demonstrated 
coupling to , which is equivalent to energy in the rotat- 
ing wave approximation. However, it remains a challenge 
to satisfy all of the requirements for achieving a practical 
QND measurement of individual quantum jumps. 

As a result, it would be desirable to detect energy 
quantization in a mechanical oscillator using only the 
presently-existing resource of a detector that couples 
weakly to energy. This is the goal of this Letter. Similar 
to previous studies, we consider the QND measurement 
of the energy fluctuations of a dissipative mechanical res- 
onator, but now consider the case where the mechanical 
resonator is strongly driven. This drive will result in a 
large average number of quanta in the resonator, n 3> 1. 
Our focus will be on the fluctuations of energy about 
this average value; in the low-temperature limit, these 
fluctuations are completely quantum in nature and re- 



flect the discreteness of the oscillator's energy. As the 
magnitude of this "phonon shot noise" scales with the 
magnitude of the applied drive, one can make it large 
enough to detect even if the detector-oscillator coupling 
is too weak to resolve individual Fock states. We show 
that this is possible even in the presence of strong cav- 
ity cooling, which is necessary to ensure that the phonon 
shot noise dominates the thermal noise. We also analyze 
the fundamental backaction of the measurement, which 
manifests itself as mechanical frequency noise that can 
obscure the intrinsic quantum energy fluctuations. We 
show that if the mechanical resonator is driven on reso- 
nance, the measurement is backaction evading, having a 
formal equivalence to a single-quadrature position mea- 
surement [5]. As such, there is no fundamental quantum 
limit on the continuous monitoring of phonon shot noise. 

Our study also reveals new physics associated with 
the energy fluctuations of a driven quantum resonator, 
namely that higher moments of the fluctuations arc far 
more sensitive to the difference between the quantum 
and classical limit than the second moment. The sec- 
ond moment in the zero-temperature quantum regime 
has a form that follows directly from the corresponding 
high-temperature, classical result: one simply takes the 
classical expression and replaces the temperature T by 
hujM /'2kB, where ujm is the mechanical mode frequency. 
We show that this simple correspondence does not hold 
for higher moments of the energy fluctuations. 

Model. — We consider a generic setup in which a de- 
tector is weakly coupled to the energy of a damped, 
driven mechanical resonator. We let F denote the de- 
tector quantity which directly couples to the oscillator 
number operator fi — c^c, and / denotes the detector 
quantity that is directly monitored. We take h = 1. The 
Hamiltonian has the form: 

H = wmc^c + H^-f (e"^"*c -t- h.c.) + H^ot + F ■ n (1) 
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Here describes the damping (at a rate 7) of the me- 
chanical resonator by a thermal oscillator bath, / is the 
magnitude of the coherent oscillator driving force (fre- 
quency 100= ujm + and Hdct is the detector Hamilto- 
nian. In the absence of thermal noise, the coherent oscil- 
lator drive would yield an average number of mechanical 
quanta (n) = n = 4/^/(4^^ +7^)- A concrete exam- 
ple of such a system is the "membrane-in-the-middle" 
optomechanical system discussed in Refs. [11[S1|7]- The 
detector here is a driven optical cavity whose frequency 
uju depends quadratically on the displacement of a di- 
electric membrane placed in the cavity; this membrane is 
the mechanical oscillator. Within a rotating-wave ap- 
proximation, one obtains the desired coupling to the 
oscillator's energy, with the input operator given by 
F = Aficn-v, where Ucav is the cavity photon number 
and A = {(PiUfi/dx'^)xlp^ is the quadratic optomechanical 
coupling. The output operator / here would correspond 
to the error signal in a Pound-Drever-Hall scheme used 
to monitor the cavity frequency (and hence n). 

For the weak coupling situation we focus on, linear re- 
sponse theory applies, implying {I{t)) = \{n{t)), where 
A is the detector response coefficient. We focus on the 
experimentally-relevant limit where the oscillator energy 
evolves slowly enough that we can ignore the frequency 
dependence of A. A straightforward calculation yields 
that the symmetrized noise spectrum of / at low fre- 
quencies is given by: 

Sll[io] = {SnnSmp + Snn[^] + Snn.BAi^]) (2a) 

SnnH = n{l + 2nti,)ip[uj + d]+p[iu~S]) (2b) 

where p[lu] = (7/2)/ (w'^ -1-7^/4), and 7 is the total damp- 
ing of the mechanical oscillator. We have anticipated us- 
ing a large magnitude oscillator driving force, and have 
thus only retained oscillator terms proportional to n. The 



first term in Eq. (2a I, Snn.imp, describes the intrinsic 



output noise floor of the detector (i.e. the imprecision 
noise), while the second term corresponds to the am- 
plified number fluctuations Snn of the mechanical res- 



onator. These number fluctuations (c.f. Eq. (2b)) have 



both a quantum shot-noise part which is non- vanishing at 
zero-temperature and a classical part proportional to the 
effective number of thermal quanta in the oscillator, nth- 
This thermal contribution corresponds to the amplifica- 
tion of the thermal force noise driving the oscillator by 
the coherent drive. Finally, the third term in Eq. ( pa] ) de- 
scribes a backaction contribution to the output spectrum: 
by virtue of the detector-oscillator coupling in Eq. ([T]), 
fiuctuations of F will result in fluctuations of the me- 
chanical oscillator frequency, thus enhancing its number 
fluctuations and yielding extra noise in the output. 

Resolving phonon shot noise. — We first discuss the 
resolvability of the phonon shot noise ignoring the ef- 
fects of backaction. A first requirement is to have the 
quantum shot noise contribution to Snn overwhelm the 



classical, thermal contribution. We thus require a cold 
oscillator, nth 1- In the optomechanical setup, this 
could be achieved by using a second optical mode whose 
linewidth is smaller than ujm to laser-cool the mechanical 
mode, as discussed in Refs. [SHIP]: such a simultaneous 
use of different optical modes for cooling and measure- 
ment can be achieved in the device of Ref. [?]• The use 
of laser-cooling comes at a price: it increases the total 
mechanical oscillator damping compared to its intrinsic 
value 7o. This in turn reduces the resolvability of the 
oscillator peak in the output spectrum. One has the sim- 
ple relation 7 = 7o("-bath/"-th), where ribath is the bath 
temperature (expressed as a number of quanta) . 

Resolving the peak associated with 5„n[a;] in the out- 
put spectrum corresponds to continuously measuring the 
phonon shot noise. Even without backaction effects, it is 
not clear that this can be accomplished, given both that 
strong cooling is required, and that the intrinsic detector- 
oscillator coupling is weak. Consider the case (5 = 0, 
which gives a maximal oscillator-induced peak in Sjj [uj] . 
As a measure of the resolvability of the phonon shot noise, 
we consider the peak-to-noise ratio S = Snn[0]/ Snn. imp- 
In the interesting limit nth — > 0, fi ^ 1, we find: 

4n/7 4n 
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^ = 8nnthS(°) (3) 
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Here, S(°) = l/(27onbathS'„„,imp) is the signal-to-noise 
ratio introduced in Ref. [1] associated with resolving 
a quantum jump of the mechanical resonator from its 
ground to first excited state. The fact that S oc nth re- 
fiects the increase in 7 associated with cavity cooling. 
We see that with a suitably strong mechanical drive, 
S > 1 even though the measurement strength is weak 
(i.e. E(o) < 1) and though a significant amount of back- 
action cooling is required (i.e. nth <C 1 ^ ?^bath)- 

We can easily apply our general analysis to the 
membrane-in-the-middle device demonstrated in Ref. [7] . 
This device is similar to those described in Refs. [S], 
but achieves values of d^uJu/dx^ up to three orders of 
magnitude greater while maintaining negligible optical 
absorption [7]. The parameters of this device are listed 
in Table 1. We assume the device is pre-cooled inside a 
cryostat to a bath temperature T = 300 mK; this would 
ensure that laser cooling to nth < 1 is readily feasible 
[21 [TO] . We also assume the membrane is driven to an 
amplitude of 2 nm. This is below the onset of dynam- 
ical bistability in these devices [H], and corresponds to 
mean phonon number n = 4.76 x 10^^. The resulting 
peak-to-noise ratio is SS > 20. Thus, by combining the 
measurement scheme presented here with the enhanced 
d'^oju/dx'^ demonstrated in Ref. [7], it should be possible 
to observe energy quantization in an existing device. This 
is in sharp contrast to the proposal for detecting individ- 
ual phonons in an optomechanical device, which would 
require substantial improvements to the membrane, the 
cavity, and their coupling. 0] 
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TABLE I; Optomechanical device parameters. m,ujm,Q' membrane mass, resonance frequency, and quality factor. F,L,Pin- 
cavity finesse, length, and input power. T: bath temperature. With the exception of T and nth these parameters are those of 
the device demonstrated in Ref. |7l[ll]. They allow a peak-to-noise ratio SS = 20.6 for the phonon shot noise measurement. 

m: 40 ng 1 MHz Q: 1.2 x 10^ F: 5 x lO" L: 67 mm Pin: 10 ^M T: 300 mK 3 MHz/nm^ na-,: 0.2 



precision noise of the detector (first term in Eq. (2a)) 



Backaction. — We now consider the effects of the fun- 
damental measurement backaction. These result from 
fluctuations of the input operator F, which lead to fre- 
quency noise of the mechanical oscillator and enhanced 
mechanical energy fluctuations. One might expect that 
to prevent this backaction contribution from obscuring 
the phonon shot noise signature, there will be a limit 
to how large one can make n, and hence a limit on the 
maximum signal to noise ratio S in Eq. ([3| . A similar sit- 
uation arises in the continuous monitoring of zero-point 
position fluctuations, where the analogous peak-to-noise 
ratio 5 < 3 [H]. 

To analyze the effects of backaction in the large n limit 
of interest, we let d denote the fluctuations of the me- 
chanical oscillator annihilation operator from its average 

value: c(t) = e"*"°* (^y/fie'f^ + d(t)j . In the rotating 
frame, the mechanical oscillator Hamiltonian takes the 
form Hm = —6d^d. Further, to leading order in n, the 
oscillator number operator becomes h(t) ~ n -|- \/^X{t) 
and the oscillator-detector coupling Hamiltonian Hi^t = 
y/nPx, where X{t) = (^e-'f^d{t) + e'f'd^t)y Measuring 
the oscillator energy thus corresponds to a measurement 
of the "position" operator X , and a continuous measure- 
ment of the phonon shot noise to a measurement of the 
zero-point fluctuations of X. As with standard contin- 
uous position detection, there will in general be a back- 
action associated with this measurement, as X does not 
commute with Hm- However, in the special case of a 
resonant oscillator drive (i.e. S — 0), Hm vanishes. X is 
thus trivially a QND observable that can be measured 
without a backaction limit. Thus, for an on-resonant 
drive, and to leading order in n, there is no backaction 
associated with the phonon shot noise measurement; we 
thus anticipate that there is no backaction-induced quan- 
tum limit on S. Formally, the measurement of phonon 
shot noise in this large-n limit corresponds to a "QND in 
time" measurement of a single motional quadrature [5]. 

To verify that there is indeed no backaction limit, we 
need to describe the leading non-vanishing backaction 
contribution when 5 — 0. This arises from the term 
(vd in n that was neglected above. Treating the fluctu- 
ations of F as being Gaussian (as is appropriate for the 
weak couplings we consider), we find that the leading- 
order backaction-driven number fluctuations are given 

by Snn,BA[0] = ^ (^) • The total added noise of 
the measurement includes this contribution, plus the im- 



The total added noise can be represented as an ef- 
fective number of thermal oscillator quanta /ladd via 
"add = (5„n,imp + S'«n,BA[0]) /(2S'„„[0,nth = 0]). As- 
suming our detector has quantum-limited noise proper- 
ties and thus optimizes the Heisenberg noise inequality 



we find that for a fixed value of 



(4) 



Snn,inipSFF > 1/4 

n, the minimal possible value of nadd is given by: 

3 

min 16 (2n)^''' 

where the minimum is achieved for an optimal measure- 
ment strength satisfying Snn,iTnp ■ 7 = (32n^)^/'^. We 
thus see that even with the inclusion of backaction ef- 
fects, the added noise of the measurement (referred back 
to the oscillator) can be made arbitrarily small by driv- 
ing the oscillator on-resonance with a sufficient strength. 
Note that this is not true if one drives the oscillator off- 
resonance. In this case, even in the limit fi — > oo, the 
added noise riadd cannot be made smaller than |(5/7|. The 
device described in Table 1 is not optimally coupled, but 
backaction effects nonetheless only yield riadd — 4x 10"^. 

Finally, in the optomechanical realization of our 
scheme, one must also consider the backaction mecha- 
nism analyzed in Ref. |13j . This mechanism is absent 
in the ideal case of a one-port cavity; in the more gen- 
eral case of a two-port cavity, the constraint set by this 
mechanism on the cavity damping is much weaker for 
our scheme than for a single-phonon measurement, as our 
scheme utilizes much smaller optomechanical coupling. 

Non-classical higher moments. — While the phonon 
shot noise described by the zero-temperature limit of 
Snijci;] is a completely quantum phenomenon, its form is 
not so different from the classical, high-temperature pre- 
diction for the energy fluctuations of a driven oscillator. 
The quantum, zero-temperature limit of S'„„ [w] is simply 
obtained by taking the corresponding high-temperature 



classical expression (second term in Eq. (2b)), and re- 



placing nth with 0.5. In other words, the quantum shot 
noise is identical to the classical expression evaluated at 
a temperature T = Sa;/2fc_B. The energy fluctuations of a 
driven oscillator are however not Gaussian in either the 
classical or quantum limit; thus, a natural question is 
whether such a simple correspondence between the two 
limits also exists for the higher moments. We now show 
that this is not the case by calculating the third moment 
of the driven oscillator's energy fluctuations. 

For simplicity, we consider the fluctuations of the time- 
integrated phonon number rh = dt'h{t') in the long- 
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time limit. Experimentally, one would attempt to mea- 
sure this quantity by time-integrating the detector out- 
put I{t); its mean is simply given by (m) = {n)t. In con- 
sidering the second and higher moments of m quantum 
mechanically, care must be taken to account for the fact 
that n(t) operators at different times do not commute. 
One must treat the measurement quantum mechanically 
(not just the oscillator) to extract the proper definition 
of the higher-moments as measured in the experiment. 

The above ordering problem has been addressed exten- 
sively in the study of the counting statistics of electron 
transport through phase coherent conductors. Several 
different models of an ideal measurement yield the same 
definition for the higher moments, the so-called Keldysh 
ordering [Ml [15]. For the second central moment, one 
obtains in the long-time limit {{rh^)) = 5„„[0]t, as would 
be obtained directly from the definition of rh. However, 
for the third central moment, a far less obvious answer is 
given {5h = fi — (n) ) : 

{{nV^))K = f dh f dt2 f dhg{6h{h)5h{t2)5n{h)) 
Jo Jo Jo 

5 - ^ [1 - 9{ti - t2)eih - t2)] (5) 

The classical definition would simply have g = 1; in- 
stead, the quantum answer involves neglecting contribu- 
tions to the integral where the middle n operator ap- 
pears at the earliest time. This ordering follows from 
a consideration of how the measured operator n influ- 
ences the evolution of the density matrix of the detector 
[l4l 115] . Importantly, one can derive Eq. ^ explicitly 
for the concrete optomechanical realization of our sys- 
tem, where the mechanical mode energy modulates the 
frequency of a resonantly-driven single-port optical cav- 
ity. To see this, we assume that the output field from 
the cavity is subjected to homodyne detection, and that 
the mechanical damping is much smaller than the cavity 
damping. Using standard input-output theory [12] . the 
output field from the mixer used in the homodyne setup 
will have the form b{t) — /3 + Bn{t), where /3 parameter- 
izes the large magnitude of the classical reference beam 
used, i? is a constant proportional to the optomechanical 
coupling, and we omit terms describing shot noise. The 
intensity / = b^b is then measured using a photodetector; 
SI = wb — is the output of the measurement. 

Experimentally, one would try to extract the third 
moment of J dtSI{t) from the observed third moment 
of the number of photons detected in the output port 
of the homodyne setup within a given time-interval. 
To calculate this, one has to use Glauber photodetec- 
tion theory[16l. In particular, the third-order inten- 
sity correlation will be given by (: I{ti)I{t2)I{ty,) :) — 

{h\ta)bHh)bHt,)b{tS{tb)h{ta)) , where < h < t, de- 
notes the time-ordered listing of ti,t2 and t^. By thus 
calculating the leading contribution (in /?) to the third 




FIG. 1: (a) Third moment {(m^)) of integrated energy fluctu- 
ations resulting from phonon shot noise, as a function of the 
drive detuning 5. (b) Zoom of plot in (a), showing that at 
low temperatures, {{m"^)) can become negative in the quan- 
tum limit. 

moment of the number of detected photons in the out- 
put port of the homodyne setup, one directly recovers the 
Keldysh ordering described by Eq. ([5|. To our knowl- 
edge, this is the first demonstration of how the Keldysh 
ordering arises in a realistic measurement setup. 
Applying this definition to our system, we find: 



// 3\\ 6ni 



3 ~ 4(6/ j)^ + 16nth(nth + 1) 



(6) 



We have neglected a purely thermal contribution to 
Eq. (|6| which is independent of n. In the high- 
temperature classical limit (i.e. nth 3> 1), the skewness 
of rh is always positive, and the dependence on the de- 
tuning S of the mechanical drive enters only through the 
oscillator's susceptibility. In contrast, in the quantum 
limit nth — ^ 0, the skewness can be positive or negative 
depending on S. Thus, the third moment of energy fluc- 
tuations has a much greater sensitivity to whether one is 
in the quantum or classical limit: the quantum expression 
does not simply correspond to evaluating the classical ex- 
pression at an effective temperature T = Hlom /'^kB- Note 
that the use of the proper Keldysh ordering is crucial to 
obtaining Eq. If one simply took 5 = 1 in Eq. ([5|, 
one would get the correct high-temperature limit, but 
would erroneously find that the zero-temperature expres- 
sion corresponds to the classical expression evaluated at 
T — huJAi/^ks- Note that the full distribution of the 
driven oscillator's energy fluctuations can be directly ob- 
tained using the approach of Ref. T7] . 
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